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ABSTRACT 



Context. Hard X-ray (HXR) spikes refer to fine time structures on timescales of seconds to milliseconds in high-energy HXR emission 
profiles during solar flare eruptions. 

Aims. We present a preliminary statistical investigation of temporal and spectral properties of HXR spikes. 

Methods. Using a three-sigma spike selection rule, we detected 184 spikes in 94 out of 322 flares with significant counts at given 
photon energies, which were detected from demodulated HXR light curves obtained by the Reuven Ramaty High Energy Solar 
Spectroscopic Imager (RHESSI). About one fifth of these spikes are also detected at photon energies higher than 100 keV. 
Results. The statistical properties of the spikes are as follows. (1) HXR spikes are produced in both impulsive flares and long-duration 
flares with nearly the same occurrence rates. Ninety percent of the spikes occur during the rise phase of the flares, and about 70% occur 
around the peak times of the flares. (2) The time durations of the spikes vary from 0.2 to 2 s, with the mean being 1.0 s, which is not 
dependent on photon energies. The spikes exhibit symmetric time profiles with no significant difference between rise and decay times. 
(3) Among the most energetic spikes, nearly all of them have harder count spectra than their underlying slow-varying components. 
There is also a weak indication that spikes exhibiting time lags in high-energy emissions tend to have harder spectra than spikes with 
time lags in low-energy emissions. 
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1. Introduction 

Solar flare emission at sub-second timescales was reported in 
hard X-ray (HXR) observations by satellite-borne HXR spec- 
trome ters i n the 70s a nd 80s (jy an Beek, de Feiter & de Jagei , 
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1989 *). Examining HXR flares observed by the Solar Maximum 
Mission (SMM ) with time resolutions of 128 ms and 10 ms, 
iKiplinger etaP (Il983h found that 53 out of nearly 3000 flares 
produce several hundred fast spikes with durations as short 
as 45 ms. These energetic flare bursts on short timescales are 
believed to be nonthermal in nature, and their temporal and 
spectral properties place constraints on the physical nature of 
the source. 

Several physical mechanisms are proposed to produce HXR 
spikes as rapid and short-lasting enhancement over the slow- 
varying underlying emission. Magnetic reconnection is believed 
to provide the means of converting energy stored in the mag- 
netic fields into thermal and kinetic energies (e.g^, CarmichaeL 
1964; Sturrock, 1966; Hiravama, 1974; Kom & Pneumanlil97a 
Priest & ForbesL i2000: Aschwanden . .2004.) . Particles may be 



accelerated to high energies at the reconnection site in the 
corona and then deposited at the chromosphere emitting HXR 
radiation by bremsstrahlung. It is known that along the length 
of the pre-connection current sheet, a tearing mode insta- 
bility may occur to trigger reconnection and form magnetic 



islands dFurfli, KiUeen, & Rosenblufll 119631; ISturrocki Il966l) . 
Formation of the magnetic islands by tearing instability and 
subsequent interactions between the se magnetic i slands, namely 
the dynamic magnetic reconnection (iKliem et al.i, 2000, and ref- 
erences therein), are cons idered to account fo r fast variations 
in nonthermal emissions (" AschwandenL l2002l and references 
therein). Kaufmann (1996) also proposed that in the magneti- 
cally complex solar active regions, one would expect to have 
multiple, primeval explosive compact synchrotron sources flash- 
ing at dilferent times, building up what is usually described 
as the onset of the impulsive phase of the bursts. Different 
spikes may indeed originate at various s ites, as suggested by 
one mm-wave observation ( Correia et al.L [l995h . The primeval 
compact sources might originate from a number of plasma in- 
stabilities, suc h as in twisted magnetic fields, and tnagneti c 
flux networks (ISturrock & UchidaL fl98l ISturrock et all [l98l . 
Another proposed mechanism to produce fine structures in 
HXR and microwaves concerns nonthermal electron injections. 
Electron beams become unstable when fast particles of a non- 
thermal electron distribution outpace the slower ones from the 
thermal or suprathermal tail during their propagation along mag- 
netic fields lines. These unstable electron beams cause discrete 
electron injections and finally contribute to spiky HXR emis- 
sions. 

The kinematics of energized electrons and the energy de- 
pendence of their collisional interactions in high-density plas- 
mas in solar flares can be probed by accurate time de- 
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lay measurements between HXR emission at different ener- 
gies. There is a variety of the time delay mechanisms op- 
erating in solar flares that can be used as a diagnostics of 
physical parameters. In the thick-targ et model (Brwon, 1971; 
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Fisher, Canfield & McClvmont), '198 51: lEmslie & Nagal ■1985i: 
MacNeice et al., 1984; Mariska, Ems lie, & Lil l 19891) for HXR 
emission in solar flares, electron acceleration is assumed to occur 
in flaring loops at coronal heights, while HXR bremsstrahlung 
emission is produced in the chromosphere. Under this as- 
sumption, the velocity spectrum of the accelerated electrons 
causes time-of-flight differences that are expected to result 
in a delay of the lower energy HX R emission w i th re- 
spect to that in higher energies ( Aschwanden et all Il995t 
iBrown. Conway, & Aschwandenlll998l) . A different case is that 
the HXR emission in higher e nergies lags behind the emission 
in lo wer energies ( Aschwanden et all 1 1997b lOiu. Lee. & GarvL 
12004 ). Two basic models have bee n proposed for the lat- 
ter: the trap-plus-precip i tation model (iMelrose & Browiil 1976 ; 
Vilmer. Kane & TrottetL Il982t iBespalov. Zaitsev. & Stepanov . 



1987) and the second-step acceleration model ("Bai & Ramatv 



19791) . In the trap model, the lag of the high-energy emission oc- 
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curs because the collisional timescale increases with the particle 
energy. This delays the escape from the trap and the subsequent 
precipitation and, thus, the resulting thick-target HXR emission. 
In the second-step acceleration model, a second-stage process 
is invoked to accelerate super-thermal electrons to higher ener- 
gies. Both of these models produce an energy delay with a sign 
opposite to that caused by the time-of-flight effect. In particu- 
lar, whether the low-energy emission lags or the high-energy 
emission lags depends on the relative importance of these two 
mechanisms. Involved in different mechanisms during the flare 
eruptions, we study the meaningful energy-dependent time delay 
that occurs in the spikes with high temporal and spectral resolv- 
ing observations. 

The Reuven Ramaty H igh Energy Solar Spectroscopic 
Imager (RHESSI, iLin et all .2002) has observed several tens of 
thousands of flares with unprecedented temporal and spectral 
resolution since its launch in early 2002. It uses a set of nine 
rotating modulation collimators, each consisting of a high spec- 
tral resolution germanium detector. As the spacecraft rotates, the 
grids transmit a rapidly time-modulated fraction of the incident 
flux. To identify fine time structures in flare HXR emissions ob- 
served by RHESSI, a demodulation algorithm was recently de- 
veloped to remove the modulation pattern from the HXR light 
curves. The demodulation code was applied to some flares that 
were found to produce fast-varying spike s evident in HXR emis- 
sions of over 100 keV (iOiu et al. Ll20Tl hereafter. Paper I). Ten 
spikes analyzed in Paper I exhibit a sharp rise and decay with 
a duration of less than 1 s, with the designated instrument res- 
olution of 125 ms. Energy-dependent time lags are found in 
some spikes, which are consistent with the time-of-flight of pre- 
cipitating electrons estimated from imaging observations. All 
these spikes show a harder spectrum than the underlying compo- 
nents, indicating that the spectral properties are closely related to 
timescales. 

These ten spikes analyzed in Paper I are produced in a few 
randomly selected flares. It is important to conduct a systematic 
search for fast spikes using existing databases. In the present 
study, we examine all flares observed by RHESSI in 2002, and 
apply the demodulation code to flares with more than 100 data 
counts in 12-25 keV to identify HXR spikes in a few pho- 
ton energy bands. We then conduct a statistical study of spike- 
productive flares and the temporal and spectral properties of 
HXR spikes. The paper is organized as follows. Observations 
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Fig. 1. Demodulated RHESSI HXR light curves (dark) in com- 
parison with summed raw data counts (gray). The light curves 
are acquired at 25-100 keV. The demodulated light curves cor- 
respond to the count rate y-axis and the summed raw data are 
related to the count y-axis. The top panel shows an event on 
2002 March 17 that exhibits fast- varying spikes. The bottom 
panel shows an event on 2002 February 20 without significant 
fast-varying spikes. The scales of the raw and demodulated light 
curves are indicated in the left and right y-axis, respectively. 



and data analysis are given in §2. Properties of HXR spikes are 
presented in §3, followed by discussions and conclusions in §4. 



2. Observations and data analysis 

2.1. Selection of HXR spikes 

Using the demodulation algorithm as described in Paper I, we 
examined several RHESSI observations of flare bursts for HXR 
spikes, and analyzed their temporal and spectral properties. Our 
database includes all RHESSI HXR flares in the year 2002. 
According to the RHESSI flare list, 5302 flares were observed 
in 2002. Among them, we selected relatively strong events with 
peak count rates in 12-25 keV greater than 100 counts s"' and 
significant emission in the photon energy range greater than 50 
keV. This information was directly acquired from the RHESSI 
flare list. Only 322 flares satisfied these two pre-analysis se- 
lection criteria and were analyzed with the demodulation code. 
These are our final sample for the statistical study. 

We applied the demodulation algorithm to obtain high tem- 
poral resolution light curves with a designated cadence of 125 
ms. Figure 1 shows examples of high-candence flare HXR light 
curves at 25-100 keV that are demodulated from the raw data. 
A large number of fluctuations are present in the raw data set, 
most of which, however, were removed in the demodulated light 
curves. Clearly, not all flares produce sufficiently significant 
HXR spikes to be picked out by the selection rule. Figure la 
shows an event with evident fast-varying structures, the spike at 
19:27:43 UT, while for the event in Figure lb no outstanding 
fast varying spikes can be recognized in the demodulated light 
curves. One cannot tell whether a flare produces fast- varying 
spikes from undemodulated raw data. 

To identify spikes, we subtracted from demodulated light 
curves /(f) a slow-varying component Isiow(t), which is obtained 
by box smoothing /(f) over a smoothing window vv,,,,,. Wsmt is 
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Fig. 2. Spike detection algorithm using the three methods (see 
text) for the event on 2002 August 3. Left panels show the de- 
modulated light curves at 25-100 keV, the slow-varying compo- 
nent (smooth) determined by the three methods, and the residu- 
als (bottom of each panel). The smoothed lines (bottom of each 
panel) show the sigma levels using the three methods. Right 
panels show the times (x-axis) and energies (y-axis) of signif- 
icant residuals recognized by the algorithms. The numbers 1 to 
7 along the y-axis indicate energy bands of 10-15, 15-25, 25- 
40, 40-60, 60-100, 100-300, and 25-100 keV, respectively. The 
vertical line indicates the most pronounced spikes of the dia- 
grams. 



taken to be between 4 and 8 s, or 32 to 64 bins. We then ob- 
tained the residual intensity as /,(f) - I{t) - Isiow(f)- We defined 
a spike when the residual intensity, /,(?), in the 25-100 keV band 
was > Usigcr for three consecutive time bins, where cr is the stan- 
dard deviation of the residual intensity profile, and Wj/g = 3, 4, 5. 
We used two ways to derive cr. First, we derived a fixed value of 
cr independent of flare evolution, namely, cr = ctq, where ctq is 
the standard deviation of /^(f) for the entire duration of the flare, 
which is usually 2-3 minutes in our data (Figure 2a). Second, to 
minimize the effect of statistical noises associated with instanta- 
neous photon counts received by the detector during the evolu- 
tion of the flare, we also used a time-dependent cr(f) derived in a 
running box Wsig (Figure 2b-c). Derivation of Ishwit), hit) and cr 
and selection of HXR spikes by these algorithms is illustrated in 
Figure 2 using the event on 2002 August 3 as an example. The 
left panels in the figure show the demodulated light curve at 25- 
100 keV, superposed with the slow component IsIok derived with 
the three methods. The residuals /,• and cr levels (smooth line) 
are also plotted at the bottom of each panel. The right panels il- 
lustrate the times (x-axis) and energy bands (y-axis: see figure 
caption) whose residuals exceeded 3cr level. Evidently, the dif- 
ferent choices of cr by different combinations of w„„/ and Wj/g 
result in different populations of spike-productive flares, which 
will be discussed below. We remark that the selection rules likely 
give a conservative lower limit of the spike population because 
of the requirement of three consecutive time bins. 

We applied the selection rules to the 322 samples. Table 1 
gives the number of spikes and the number of spike-productive 
flares using different combinations of Wj,,,, and Wj/g, including 
the fixed erg criterion, and with n „g - 3. Depending on how cr is 
defined, the selection rule of 3cr and three consecutive time bins 
yields 40 to 90 spike -produtive flares, each producing one or two 
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Fig. 3. Demodulated RHESSI HXR light curves at varying ener- 
gies for three flares. The vertical lines indicate the spike times. 
The y-axis gives the units of the 25-100 keV light curves, and 
light curves at other energies are scaled arbitrarily. 
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Fig. 4. Number of detected spikes varies as cr level increases. 
The asterisks and the line represent observational and Poisson 
distributions, respectively. 



spikes, out of the total of over 5302 flares observed by RHESSI 
in 2002. This productivity is in general comparable with the re- 
sult by Kiplinger et al. ( 1983), who, also by applying the 3cr and 
consecutive bin criterion, detected 53 spike-productive flares out 
of 3000 flares observed by SMM. That is to say, fast-varying 
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Table 1. Spike/flare number variations with cr defined by differ- 
ent combinations of Wsmt and Wsig and with n„g = 3. 





w.™, = 4s 


Ws„, = 6s 


w,„, = 8s 


3o-o 


184/94 






Wj;g=60s 


70/57 


91/74 


103/80 


w„g=45s 


67/54 


76/60 


91/72 


Wj,j=30s 


47/39 


57/46 


65/51 


w„g=15s 


28/25 


27/23 


31/26 


Table 2. Spike/flare number variations with different Usig. 
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spikes are detected in 10-20% of flares, though observed by dif- 
ferent instruments with different methods. 

Figure 3 shows the demodulated light curves of three flares 
at varying energies from 25 to 300 keV during the flare erup- 
tions. The designated time bin is 0.125 s. These events all exhibit 
fast- varying structures on timescale of 1 s or less at most of the 
energy bins above 25 keV. Signals up to 100 keV are significant 
and are likely real features. Some of these spikes are also visi- 
ble in 100-300 keV, though less significant because of reduced 
counts level. 

Using the fixed 3(To yields 184 spikes in 94 flares. At 5(To, 
there are 54 spikes in 40 flares. At face value, these numbers 
are significantly greater than the number of spikes that could be 
produced by statistical noises if the noise distribution followed 
a simple Poisson (Figure 4). In reality, the noise distribution is 
poorly defined, and the noise dependence on instantaneous pho- 
ton counts is likely much stronger than a Poisson. Therefore, we 
performed a negative test to identify "negative" spikes by setting 
Hsig = -3, -4, -5 and using the same rule of three consecutive 
bins. The results of the negative test for a few representative se- 
lections are fisted in Table 2. Evidently, the negative occurrence 
is nearly one half of the positive occurrence for the fixed ctq. 
Even though one cannot ascribe the negative occurrence entirely 
to the statistical negative because of the way the residual signals 
are produced, it is conservative to say that about one half of the 
positive occurrence probably reflects real spike signals. 

When the time-varying cr is employed, understandably, the 
number of spikes steadily decreases with a smaller Wsig or vv,m, 
(Table 1). Table 2 clearly shows that the negative occurrence 
with the time- varying cr drops to only 10-20% of the positive 
occurrence. Therefore, most spikes picked out with the time- 
varying cr approach may be considered as real spikes free from 
statistical noises. We also note that 2/3 of the spike population 
via the time- varying cr approach are included in the population 
from the fixed Sctq selection rule. 

The selection with the time-varying cr does not only change 
the total number of spikes but would, probably, change the dis- 
tribution of these spikes during the evolution of the flare. For 
example, the fixed ctq regardless of the evolution of the flare 
would misrepresent the noises at different evolution stages, cro 
tends to be smaller than the real standard deviation determined 
by statistical counts during the flare maximum but larger than the 
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Fig. 5. Peak count rate distributions at various criteria denoted 
in the figure. The upper thick solid line indicates our 322-flare 
sample. The profile is obviously steeper than all other profiles. 



standard deviation during the rise-and-decay phases of the flare 
when the counts are lower than at the peak. Therefore, the con- 
stant cTo rule might include noisy signals at the peak of the flare, 
while overlooking real signals during the rise-and-decay phase. 
On the other hand, by using the time-varying cr, the sample is 
partly biased against high counts, or the peak phase of the flare. 

Because we do not know the true noise distribution, we 
examined and compared properties of spikes and associated 
flares selected with the following different rules: fixed cro, time- 
varying cr with Wsmt = 8 s and w„g = 60 s, and time- varying cr 
with Wsmt = 4 s and Wsig ' 30 s. The spikes were selected with 
Hsig = 3. Specifically, we illustrated the distribution of the events 
with respect to the flare magnitude, flare duration, and spike oc- 
currence time, and examine statistical properties of the spikes in 
terms of their evolution and energetics. 

2.2. Event distribution 

We compared spike-productive flares with flares that do not pro- 
duce fast-varying spikes to examine whether they are statisti- 
cally different populations in terms of flare magnitude and dura- 
tion. Figure 5 shows the peak count rate distributions of all 322 
flares in the sample with a peak count rate exceeding 1 00, and the 
spike-associated flares from four different selection rules: 5£ro, 
3cro, 3cr at Wsmt = 8 s and w^ig - 60 s, and 3cr at w.,,,,, - 4 
s and Wsig = 30 s. Figure 5 shows that the distribution for the 
entire sample of 322 flares is steeper than that for the spike- 
associated flares, suggesting that more intensive flares have in 
general a greater chance to produce spikes. There is no notable 
difl'erence in the distributions of the four spike-productive flare 
populations, suggesting that the majority of selected spikes are 
not dominated by noise that grows with photon counts. 

We also examined whether impulsive flares, compared with 
gradual events, are more productive in fast-varying spikes. 
Figure 6 displays the rise time distributions for spike-productive 
flares in comparison with that of all 322 flares. Here we define 
the flare rise time as the time difference between the flare peak 
time and the start time, both referring to emission in 12-25 keV 
provided by the RHESSI flare list. No evident difference is found 
between these distributions, suggesting that HXR spikes can oc- 
cur in both impulsive and gradual events. 

For the spike-productive flares, we also studied at which time 
the spikes occur with respect to the evolution stage of the flare 
bursts. For convenience, we defined a normalized spike occur- 
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Fig. 7. Occurrence time Tgpk distribution at various criteria as 
denoted in the figure. See §2.2 for details. 



rence time Tspk, which refers to the time difference between the 
spike peak time and the flare start time normalized to the flare 
rise time. Therefore, if T^pi^ < 1, the studied spike should occur 
in the rise phase of the flare; if T^pk > 1, the spike occurs during 
the decay phase. Figure 7 illustrates the occurrence time distri- 
butions for the spike populations from varying selection rules. 
For all populations, nearly all spikes occur before or during the 
peak of the flare. Clearly, the fixed ctq selection yields a large 
percentage, up to 70%, of spikes around the peak times (Tspk 
= 0.8-1.2) of the associated flares. With ever smaller w„g and 
Wsmt, the fraction of spikes during the flare peak is reduced as 
the cr around the flare peak time increases, while there are more 
spikes detected during the low-count rise phase. At this stage of 
research, as the noises of the demodulated flare Ught curves are 
not well understood and cannot be independently determined, 
we did not readily give a preference to a certain selection rule. 
Therefore, all these spike populations from difl'erent selection 
rules are analyzed in the following study. 

Fast-varying spikes were first selected using integrated 
counts in 25-100 keV. We then demodulated the light curves of 
these events in photon energy ranges of 25^0, 40-60, 60-100, 
and 100-300 keV, and searched for the highest energy in which 
spikes are detected. We denote this energy band by em- Figure 8 
displays the spike distribution with respect to e/,, . The majority of 
spikes discovered in 25-100 keV can still be seen in 40-60 and 
60-100 keV. Indeed, for each population, the peak of the distri- 



bution is at 60-100 keV. Nearly 20% of spikes can be detected 
as high as 100-300 keV. This result suggests that fast-varying 
spikes are most probably nonthermal in nature. 

All these statistical results indicate that fast-varying spikes 
are small-scale energetic events produced during the most en- 
ergetic stage of flares. In particular, they tend to occur during 
the rise phase of intensive flares. On the other hand, impulsive 
and gradual flares have an equal chance to produce fast- varying 
spikes. 

In different phases of flare eruption, energy release, and en- 
ergy transport, the relevant physical processes are quite differ- 
ent. Theoretically, energy release occurs preferentiaUy during 
the flare rise phase. This seems to support the result of differ- 
ent productivities of spikes in different flare phases. 



3. Properties of HXR spikes 

HXR spikes at shortest timescales are beheved to reflect single 
energy release events, such as by magnetic reconnection and par- 
ticle acceleration. We examined the energy-dependent temporal 
properties of these HXR spikes to gain insight into the nature 
of the energy release. For the detected spikes, we estimated the 
spike duration, rise and decay times, the count spectral index, 
and the energy-dependent time delays and compared them with 
large flares and previous studies. 

To obtain a statistical result of the temporal and spectral 
properties of HXR spikes, we applied a simplified analysis to 
all spikes. A more rigorous analysis was also conducted on a 
dozen prominent spikes out of the selected samples (Paper I), 
which yields results consistent with the statistical results in the 
present paper. 



3.1. Duration, rise, and decay times 

The energy-dependent duration of a spike, T(e), is estimated to 
be the time difference between the e-slope start and end times 
of the spike. We determined the rise and decay times when the 
spike net intensity, i.e., the intensity with the interpolated under- 
lying component subtracted, drops to e ' of the peak intensity 
in both the rise and decay phases. As displayed in Figure 9, T(e) 
ranges from 0.25 to 2.0 s, or 2 to 16 time bins, with an average 
being 0.98±0.47 s in 25-100 keV for all HXR spikes. The same 
analysis was applied to light curves in different energy bands, 
resulting in a mean duration of 1.08+0.39 s, without significant 
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ples of HXR count spectra of integrated spike flux during the 
spike time in comparison with the count rates spectra of un- 
derlying components. Solid lines show the least-squares fit of 
the spectra to a power-law distribution, and the exponents and 
their uncertainties from the fit are noted in the figures. The y- 
axis gives the units of the spike intensities, and the intensities 
of the underlying component are divided by a number for the 
clarity of display. SI is the example we discussed in Figure 2. 



difference in different energies. Therefore, our studies show that 
the spike duration is nearly independent of photon energies. 

HXR spikes appear to have quite a symmetric sharp rise and 
decay. The mean e-slope rise time of all spikes is 0.52+0.26 s, 
and the mean decay time is 0.50+0.25 s. The result indicates 
no significant difference between the rise and decay times. This 
is somewhat different from the general scenario of flares, i.e., 
HXR emission rapidly increases in the impulsive phase, and then 
decreases relatively gradually in the decay phase. 

3.2. Spectral index 

To study the spike energetics, we fit the spike count spectrum 
at >15 keV to a power-law distribution ~ e^^ to obtain the 
power-law index F, as indicated by the straight lines in Figure 
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Fig. 11. Scatter plot of spectral indices of the spikes versus that 
of the underlying components using data from 15 keV. 



10. Only data in photon energies greater than 15 keV were used 
to avoid the contamination of thermal bremsstrahlung in lower 
energies. Among these spikes, reliable least-squares fitting re- 
sults can be made for only a fraction of spikes. The spectral fit- 
ting was applied to the net spike counts integrated over its dura- 
tion. The net counts of the spike were derived by subtracting the 
underlying component, which is a mean between the pre-spike 
and post-spike counts, also integrated over the duration of the 
spike. As a comparison, we also fit the count spectrum of the 
underlying component. Figure 10 gives two examples of inte- 
grated spike flux with respect to photon energy, i.e., the spike 
counts spectrum. Also plotted are the directly measured count 
rate spectra of the underlying components adjacent to the spikes. 
Figure 10 shows the comparison of spectral indices between the 
spikes and underlying components. Note that F is the index of 
uncalibrated count spectrum, which is somewhat different from 
the index of HXR photon spectrum or nonthermal electron spec- 
trum. However, it is still a meaningful parameter when compar- 
ing the energetics of the spike relative to the underlying compo- 
nent. Figure 1 1 shows that nearly all spikes have a lower spectral 
index than their underlying components by more than a err, the 
gross error in the fitting. Because we are aware of possible ther- 
mal contribution to photons below 25 keV, we also fit the data 
from 25 keV for comparison. The fitting result is nearly the same 
as that including the 15-25 keV counts. This statistical result is 
consistent with the more sophisticated case study in Paper I, in- 
dicating that most of the spikes are more energetic, with a harder 
nonthermal spectrum by about 0.5 in count spectrum index than 
the underlying components 

It is known that flare HXR emission exhibit s a harder spec- 
trum at emission peaks than at valleys (e.g., iKiplinger et al.L 
Il983h . The spectral analysis for spikes agrees with this gen- 
eral scenario. The different spectral characteristics between fast- 
varying HXR spikes and underlying components are likely a 
result of specific physical mechanisms to generate nonthermal 
emissions in such a way that the spectral properties are closely 
related to the burst timescales. 



3.3. Energy-dependent time delay 

We analyzed the time-delay properties of all the spikes detected 
by different criteria. Four energy bands (25-40, 40-60, 60-100, 
and 100-300 keV) were used in our analysis. For each spike, 
we found the peak times T(e) at different energies relative to the 
peak time in 25-100 keV. We then performed a least-squares 
linear-fit between the mean photon energy and r(e) in logarith- 
mic scale for each spike. If the fitting error was small, we as- 
sumed that the spike exhibits a systematic energy-dependent lag 
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Fig. 12. Distributions of spike peak time lags between 60-100 
keV and 25^0 keV for spikes selected with various criteria. 
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Fig. 13. Scatter plot of the spectral index versus the time de- 
lay parameter a of spikes selected with various criteria. In each 
panel, positive a values refer to high-energy delayed events 
while negative ones refer to low-energy delayed spikes. These 
two populations are separated by the vertical dashed line, and 
the number of the spikes, and the mean count spectral index F 
and its standard deviation of each population are also marked in 
each panel. 



in its peak times. Indeed, less than 40% of the spikes evolve with 
regular time delay patterns, and both low-energy delay and high- 
energy delay events are detected (note that in Paper I, we only 
found low-energy delay events). In Figure 12, we iUusttate the 
distribution of time lags between 60-1 00 keV and 25^0 keV for 
these events. In the figure, positive time lags indicate that high- 
energy emission lags behind low-energy emission, and negative 
time lags indicate that low-energy emission lags behind high- 
energy emission. Evidently, the majority of events exhibit time 
lags shorter than 0.5 s. The mean time lag is about 0.80 s and 
-0.74 s for high-energy delays and low-energy delays, respec- 
tively. 

We additionally investigated energy-dependent time lags 
with respect to the spectral properties for a small fraction of 
spikes that exhibit a systematic time lag pattern as well as re- 
Uable spectral fitting results (§3.2). For these spikes, we derived 
spike peak times T(e) at different energies relative to a refer- 
ence energy (25-100 keV), and fit T(e) to an arbitrary function 
T(e) ~ e". A positive a indicates that peak emission in higher 
energies lags behind the peak emission in lower energies, and 
a negative a indicates lower energy emission lagging behind 
higher energy emission. The time delay parameter a for difi'er- 
ent populations of the spikes is shown in Figure 13. We divided 
the spikes into two types: low-energy delayed (negative a) and 
high-energy delayed (positive a). Accordingly, the left part of 
each panel in Figure 13 displays the spike spectral index distri- 
bution for low-energy delayed events, while the right part shows 
the distribution for high-energy delayed spikes. In the last panel, 
the sample number is two, and both of them are low-energy de- 
layed spikes. These figures show that, among these events, about 
or more than 2/3 are low-energy delayed, and the rest are high- 
energy delayed. Furthermore, it is found that, on average, high- 
energy delayed events have a harder count spectrum than low- 
energy delayed events. 

The energy-dependent time lag patterns associated with 
spike energetics may result from the acceleration and/or trans- 
port mechanisms. For example, depending on the initial pitch- 
angle distribution, when a significant amount of electrons are 
trapped in the corona, Coulomb colUsions in the trap will eject 
these electrons from the trap to precipitate and lose their energy 
in the chromosphere. The rate of Coulomb collision is energy- 
dependent, causing low-energy electrons to precipitate first, thus 
producing high-energy delayed events. On the other hand, when 



most electrons are directly precipitated in the chromosphere, the 
energy-dependent time delays are associated with the electron 
time of flight, i.e., higher energy electrons with a higher veloc- 
ity will precipitate ahead of low-energy electrons, accounting 
for the low-energy delayed events. Our preliminary result show- 
ing high-energy delayed events with a harder spectrum indicates 
that the energetics of HXR spikes are probably coupled with the 
transport eff'ects, which are, in the first place, governed by accel- 
eration mechanisms determining the electron pitch-angle distri- 
bution. A better understanding of the mechanisms that generate 
these spikes should be obtained from the knowledge of their spa- 
tial and magnetic properties. 

4. Discussions and conclusions 

We have searched all flares observed by RHESSI in 2002 for 
fast-varying HXR spikes, finding that at least 20% of the flare 
bursts with count rates exceeding 100 counts s"' in 25-100 keV 
produce HXR spikes, mostly during the rise phase of the flares. 
These spikes have timescales of 0.2-2 s in photon energy range 
of up to 300 keV. The main results are as follows. (1) Both im- 
pulsive and long duration flares can produce HXR spikes with 
nearly equal production rates. Flares with high peak count rates 
are more productive in HXR spikes. (2) Almost all spikes occur 
in the rise phase of the flares, and a large percentage, up to 70%, 
of spikes are produced at or about the flare peak times. (3) The 
mean duration of spikes is about 0.9-1.0 s, independent of pho- 
ton energies. The rise and decay times of spikes are shown to be 
almost the same. This differs from ordinary flares that usually 
have a longer decay phase dominated by thermal emission. (4) 
Most of the spikes can be detected in very high energy bands up 
to 100-300 keV. The HXR spectra of spikes are harder than those 
of the underlying slow-varying components. This fact implies 
the nonthermal origin of spikes. (5) Evident energy-dependent 
time lags are present in a fraction of spikes, indicative of time-of- 
flight or Coulomb collision efl'ects. It is also shown that, on aver- 
age, spikes lagging in high-energy emissions have harder spectra 
than spikes exhibiting lags in low-energy emissions. 

To understand the mechanisms for HXR spikes reported in 
this study, more observational investigations are needed with 
substantially improved capabilities or techniques to conduct pre- 
cise spectral and imaging analysis at the timescales of the spikes. 
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